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Abstract 

In supersymmetric scenarios where the scalar tau lepton is stable or long-lived, 
a search for a decay mode x*^ ~^ "^A* a-t the LHC has a good sensitivity to the flavor 
mixing in the scalar lepton sector. We demonstrate that the sensitivities to the 
mixing angle at the level of sin 623 = 0.15 are possible with an integrated luminosity 
of 100 fb^^ if the total production cross section of supersymmetric particles is of 
the order of 1 pb. The sensitivity to the mixing parameter can be better than the 
experimental bound from the r — > /i7 decay depending on model parameters. 



If new physics contains a charged stable particle, such as the scalar tau lepton (f) in 
super symmetric (SUSY) models, it provides a very clean signal at the LHC experiments. Once 
they are produced, most of them penetrate detectors and leave charged tracks just like muons. 
By measuring the velocity at the muon system, we can easily distinguish from the muon 
background. A very precise mass measurement is possible by combining with the momentum 
measurements 

Scenarios with such a charged stable or long-lived particle have sounded exotic and regarded 
as alternative possibilities. However, recent serious considerations of SUSY phenomenology 
have shown that it is indeed theoretically motivated Q. The presence of such particles does not 
immediately contradict with cosmological history. There are interesting cosmological scenarios 
and even motivations for such a long-lived particle |5|, ^. If it is the case, we will have new 
kinds of signals in new physics search experiments. 

In this paper, we propose a search strategy for flavor mixing in the scalar lepton sector in 
the stable (or long-lived) f scenario at the LHC. In the presence of the flavor mixing, we will 
have a decay mode of the neutralinos ~^ '^l^- By looking for sharp peaks in the f-^ invariant 
mass, we show that we will be able to discover lepton flavor violation for r(xi ffj,)/T{xi ~^ 
tt) ~ tan^ 023 ^ 10~^, where 623 is the slepton mixing angle. 

There have been many studies on lepton flavor violation at the LHC and e~^e~ colliders 
assuming the neutralino to be the lightest SUSY particle. The possibility of observing e^fi^ 
+ missing Et final states at e~^e~ colliders has been pointed out in Ref. The correct 
treatment of the process including quantum interference (slepton oscillation) has been studied 
in Ref. and discussion has been expanded to the LHC experiments and CP violation in 
Ref. 1^]. Following those papers, LHC studies on searches for decay processes X2 ~^ ^t^JXi 
with i ^ j have been done in Refs. |1C]. The sensitivities of 0(0. 1) for mixing angles have been 
derived in various SUSY models. 

Lepton flavor violation in the long-lived f scenario has also been studied. In Ref. |11], the 
decay of f into e or /x and a gravitino is studied under an assumption that a significant number 
of f's will be collected at the LHC or future linear collider experiments by placing a massive 
stopper material close to the detectors ]12U. A linear collider study with long-lived f has also 



been done in Ref. [13| where it is proposed to search for lepton flavor violating flnal states such 
as (e'^T^f^)f^ through slepton pair production processes. Very good sensitivities as well as 
s'mO ~ (a few) X 10^^ are reported in both of the works. We study in the following the LHC 
signals of slepton flavor mixing without new detectors or future colliders. Therefore it serves 
as the first search strategy that can be done immediately after the LHC starts if f is stable or 
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long-lived. 

In order to estimate the sensitivity, we performed the following Monte Carlo simulation. We 
used a model of Ref. ^ where the spectrum of the SUSY particles are parametrized by four 
quantities: 

-^gaugino(= -^^/^s)) -^messi -^mess- (1) 

The fi parameter and Mgaugino control the Higgsino mass and the gaugino masses, respectively. 
The messenger scale Mmess and the number of messenger particle A'mess determines masses of 
scalar particles relative to the gaugino masses. We have chosen two parameter points where f 
is the lightest SUSY particle (except for the gravitino): 

Model I : /u = 300 GeV, Mgaugino = 900 GeV, M^css = 10^° GeV, iY^css = 1 , (2) 
Model II : n = 500 GeV, Mgaugino = 900 GeV, M^ess = 10^ GeV, iVmess = 1 • (3) 



The f and neutralino masses in Model I are mf = 116 GeV and (m^o, m^o, m^o, m^o) = 
(187 GeV, 276 GeV, 306 GeV, 404 GeV). In Model II, mf = 157 GeV and (m^o, m^o, m^o, 
m^o) = (194 GeV, 346 GeV, 505 GeV, 525 GeV). We chose a parameter with heavier Higgsinos 
in Model II. 

With the SUSY spectra, we have generated 40,000 SUSY events for each model by using 



the Herwig 6.50 event generator ||l^ with the CTEQ5L parton distribution function |15|. This 
corresponds to an integrated luminosity of 33 fb"^ (46 fb~^) at the LHC for Model I (Model 
II). We set the mixing angle of the right-handed sleptons to be sin023 = 0.33 with which the 
branching ratio of the lightest neutralino is T{xi ~^ ^f^)/^iXi ~^ t^t) — O-l- Heavier neutralinos 
do not have significant branching ratios for the ~^ decays because the amplitudes are 
suppressed by the Yukawa coupling constant of the muon. With SUSY spectra with the lightest 
neutralino being almost the Bino and the lighter f to be almost right-handed (which is the case in 
the above two models), the following method is not sensitive to mixings in left-handed sleptons. 



The events are passed through a detector simulator AcerDET 1.0 [16| where muon momenta 



are smeared according to the resolutions of the ATLAS detector. We have also smeared the 



momenta and velocities of f's according to the resolution obtained in Ref. |17]; 



-^ill = kip®k2\l + ^®—, (4) 
p V p p 

where ki = 0.0118%, k2 = 2% and = 89%. The momentum p is in GeV. The resolution of 
the velocity is 

^ = 2.8% xp . (5) 
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Figure 1: The f-fi invariant mass distribution in Model I. The shaded histogram in the left 
panel is the estimated background, i.e., the f-e invariant mass distribution. By subtracting the 
estimated background, we obtain the histogram in the right panel. 



We have ignored the dependence of the resolutions. Also, in the following analysis, we assume 
that the f mass is known with a good accuracy by the method of Ref. 

We have followed the strategy of Ref. Q for the identification of f . We require the candidate 
tracks to be within |r/| < 2.4, Pt > 20 GeV and /37meas > 0.4. The cut on the measured 
velocity ensures f to reach the muon system. A consistency condition: — /3meas| < 0.05 is 
imposed, where /?' is a velocity calculated from the momentum = y'pV(p^~+^"|0 ) • By also 
requiring the measured velocities of at least one candidate f to be within 0.4 < (3jraeas < 2.2 for 
each event, this selection strategy reduces background from mis-identified muons to a negligible 
level We therefore ignore in the following analysis the background from the standard model 
processes as well as from muons in SUSY events. 

In order to look for lepton flavor violating neutralino decays, we selected events with only 
one isolated muon with Pt > 20 GeV and at least one opposite-sign f candidate. If there are 
two opposite-sign f-^i pairs, we use both of them for the analysis. The invariant mass Mf^ is 
calculated for each candidate event. 

The invariant mass distribution is shown in the left panel of Fig. |l[ where we can clearly 
see a peak at the lightest neutralino mass (187 GeV). The shape and normalization of the 
background distribution can be obtained from the M^e distribution directly from the data 
(shaded histogram). By subtracting those estimated background, we obtain the histogram in 
the right panel where we see that the background is successfully subtracted. Therefore we 
can reliably use the Mfe distribution as an expected background. We can also find an excess 
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around masses of the hghter Higgsino-Hke neutrahno (~ 276 GeV). Most of the background 
originates from the ~^ ''"t decay followed by r ^ ixvi?. Although the signal region is the 
kinematic endpoint of this background for each neutralino, the background f /u pairs from heavier 
neutralinos fall into the signal region. There are also backgrounds from leptonic decays of W 
bosons. Numbers of such background events depend on cascading pattern of heavy SUSY 
particles. 

We fitted two peaks in the right panel of Fig. |^ with the gaussian function and defined 
the signal region to be the la region around the peaks; |Mf^ — 185.2 GeV| < 3.8 GeV and 
|-/Vff^ — 276 GeV| < 10 GeV. In the case where there are not enough events to find the neutralino 
masses by the Mf^ distribution, one should look for edges in the invariant mass of f and r- 
jet, Mfj^, for the neutralino mass measurements as is done in Ref. There are S + B = 
584 events in the signal region whereas the number of the expected background in the signal 
region is -B = 374. Therefore we obtain 9a excess with 33 fb~^ of data.0 Normalizing to the 
integrated luminosity of 100 fb~^, the number of needed signal events for 5a discovery to be 
181, corresponding to sin 023 > 0.18. 

We repeat the same analysis for the Model II and the result is shown in Fig. ^. There is only 
one peak associated with the Bino-like neutralino because the number of Higgsinos in cascade 
decays is reduced and the branching ratios of -B(x3 4 ~^ ^t) are suppressed with the relatively 
heavy Higgsinos. A slightly better sensitivity than Model I is obtained. In the signal region, 
|Mf^ - 193.0 GeV| < 6.9 GeV, we find 5 + 5 = 539 and B = 238 for 46 fb'^ Normalizing to 
100 fb~^ of data, we obtain the 5a sensitivity to be sin 023 > 0.15. 

If no peak is found due to small mixing angles, one can put a bound on the branching fraction 
(or equivalently the mixing angle). This requires a counting of the ~^ '^'^ events that involves 
the efficiency measurement of the r identification. That will be the dominant uncertainty in 
putting the experimental bound. As far as order of magnitude is concerned, the sensitivity will 
be at the level of sin 023 ~ 0.1 — 0.2. 

The search for the r — > //7 decay has already put a stringent bound on a combination 
of various SUSY parameters involving the slepton mixings. Although a model independent 
comparison is not possible, we can get a sense of sensitivities to the mixing parameter by 
calculating the r — > /i7 branching ratio with a particular parameter set. We have done that in 

*This level of excess is somewhat optimistic given that we know the correct location of the peaks. In the actual 
experimental situation, the peak locations (the neutralino masses) will be measured by looking for the endpoint 
locations of the invariant mass Mfj^. The uncertainty of this measurement is estimated to be at most of order 
5% by taking into account the effects of fake r-jets and the uncertainties in calibration of the r-jet energies. If 
we use the central values given in Table 2 of Ref. and the 5% errors for the definition of the signal region, i.e., 
194 ± 10 GeV and 279 ± 14 GeV, we obtain about a 7a excess, where most of the significance is a contribution 
from the second peak. 
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Figure 2: The f-/i invariant mass distribution in Model II. The shaded histogram in the left 
panel is the estimated background, i.e., the f-e invariant mass distribution. By subtracting the 
estimated background, we obtain the histogram in the right panel. 
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Table 1: LHC sensitivities and comparison to the r fj.'j decay. The number of signal {S) and 
background (B) events are shown for 40,000 SUSY events and sin023 = 0.33. For this level of 
the large mixing angle, the statistical significances can be as large as lOo". The 5(T-level discovery 
with an integrated luminosity of 100 fb~^ requires the angle to be sin 023 > 0.15. The r ^ //7 
branching ratios are shown for sin023 = 0.15. 

Model I and II with sin023 = 0.15. The branching ratios are B{t ^ /i7) = 1 X 10~^ and 4 x 10"^ 
for Model I and II, respectively. Compared with the current experimental bound, 4.5 x 10^^ 
the LHC sensitivities can be much better (or worse) depending on model parameters. One should 
also note that the measurement of T{x^ ~^ t^I^)/'^{x^ ~^ t^t) the LHC will directly probe the 
slepton mixing parameter. Therefore, measuring/constraining the branching fractions of both 
processes will be important to understand the flavor structure of SUSY models. We summarize 
the results in Table |l[ 

A similar analysis will go through for a ^ search at the LHC. Also, if a linear collider is 
built in future, searches for a decay mode ~^ '^/^ through a neutralino pair production process 
may give a better sensitivity to the mixing angle as background from heavier neutralinos and 
W bosons will be under better control. 
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